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Abstract 

We study the AdS 5 x S 5 sigma-model truncated to the near-flat-space limit to 
two-loops in perturbation theory. In addition to extending previously known 
one-loop results to the full SU(2|2) 2 S-matrix we calculate the two-loop cor- 
rection to the dispersion relation and then compute the complete two-loop 
S-matrix. The result of the perturbative calculation can be compared with the 
appropriate limit of the conjectured S-matrix for the full theory and complete 
agreement is found. 
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1 Introduction 

There has been much recent progress in the effort to completely establish the AdS/CFT 
correspondence [1]. The full conjectured integrability of planar Af = 4 Super Yang- 
Mills [2] and its dual theory, the string sigma-model on an AdS 5 x S 5 target space [3], 
has been instrumental in this progress. At least for the question of gauge operators of 
infinite bare dimension, computing the spectrum has basically come down to finding a 
two-particle S-matrix [4] that can be determined for both large and small values of the 
't Hooft coupling. 

For large 't Hooft coupling, the scattering is that of string oscillators on the world- 
sheet, while for small 't Hooft coupling it more closely resembles the scattering of 
magnons on a spin-chain. Remarkably, as was shown by Beisert [5,6], the S-matrix 
is almost completely determined by the underlying superalgebra su(2|2) x su(2|2) with 
central extension, no matter what the coupling. The only part of the S-matrix that 
cannot be determined from the supergroup structure itself is an overall phase factor (the 
dressing phase), which was conjectured first in the form of an asymptotic series at strong 
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coupling [7], and then non-perturbatively [8]. First steps towards derivation of the dress- 
ing phase from Bethe ansatz were taken in the recent work [9]. The conjectured dressing 
phase makes the S-matrix crossing-symmetric [10] and passes a remarkable four-loop test 
at weak 't Hooft coupling: decoration of the Bethe equations with the conjectured phase 
modifies the anomalous dimensions starting from four loops and such a modification 
brings the Bethe-ansatz prediction for the cusp anomalous dimension [8] into agreement 
with the explicit four-loop calculation [11]. 

It is remarkable that explicit four-loop calculations in Af = 4 SYM are possible and 
it is certainly desirable to reach comparable accuracy on the string side. Currently 
state of the art is the one-loop order: quantum corrections to the energies of various 
classical string configurations have been computed in [12]. The purpose of this paper 
is to go beyond the one-loop order. Since the full AdS 5 x S 5 sigma-model [13] is quite 
complicated we make use of the simplifying limit proposed recently by Maldacena and 
Swanson [14]. 

As in [15,16], we will be interested in the world-sheet S-matrix which can be directly 
compared to the su(2|2) xsu(2|2) S-matrix [5,6] with the conjectured dressing phase [7,8]. 
The world-sheet S-matrix simplifies immensely in the Maldacena-Swanson limit, but is 
nonetheless nontrivial since the resulting sigma model is still interacting. The limit 
is taken by scaling all momenta, such that pX 1 ^ 4 is finite. The momenta of the string 
excitations then sit in the "near-flat" region, between the noninteracting BMN regime [17] 
and the classical giant magnon regime of [18]. For excitations in the near-flat region, 
although there is no sinp/2 factor in the dispersion relation like in the case for giant 
magnons, the Lorentz invariance of the BMN region is still broken by interaction terms. 
However, as we will show this breaking of Lorentz invariance is rather mild, and in fact 
can be restored if one compensates any Lorentz boost with a rescaling of the world- 
sheet coupling constant. It might be possible to argue that the S-matrix satisfies the 
usual crossing symmetry as a consequence of the usual LSZ theorems, with additional 
modifications due to this "soft" breaking of Lorentz invariance. We shall see that the 
crossing symmetry is certainly there at the level of Feynman diagrams. 

The near-flat limit also leads to a simplification of the Janik's equation [10]. The odd 
solution will still be a sum of dilogarithms, but the even phase simplifies tremendously 
and will end up being the log of a rational function of the world-sheet coupling (and so 
its contribution to the S-matrix is to multiply it by a rational function). It is simple to 
check that this function is a solution to the near-flat limit of the BHL even equation. 
The S-matrix for the various processes also turns out to be a quadratic polynomial of 
the world-sheet coupling multiplied by a common function. 

Computing quantum corrections is much simpler in the near-flat limit. The quartic 
nature of the interaction terms makes the computations similar to those found in 4 the- 
ory in two dimensions. For two-point functions, supersymmetry prevents any tadpoles 
from occuring, so there is no one-loop wave function renormalization or mass-shift. How- 
ever, at the two-loop level there are sunset diagrams which induce radiative corrections 
to the dispersion relation that agree with the predicted near-flat limit of the dispersion 
relation in [14]. 

We then consider corrections to the four point amplitudes. We will compute these 
corrections up to the two-loop level, where we will find agreement with the near-flat 
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limit of the BHL prediction. This provides the first nontrivial check that goes beyond 
the tree level AFS [19] and one-loop HL [20] dressing factor terms. In carrying out these 
computations, we will see that the final amplitudes for the different processes are very 
similar, as they must be if they are to agree with the BHL S-matrix, but the road to how 
these final amplitudes are reached can be significantly different. For example, for certain 
2^2 bosonic processes, there is a four-fermion interaction term that contributes to the 
two-loop amplitude, while in other processes this interaction term plays no role. In any 
case, the underlying supersymmetry must play a crucial part in determining the final 
structure of these amplitudes. In going from the amplitude to the S-matrix, we must 
take into account the two-loop wave-function renormalization as well as the two-loop 
mass-shift which will affect the Jacobian factor that needs to be included. 

This paper is structured as follows. In Sec. [2] we review Maldacena and Swanson's 
action for the near-flat limit. In Sec. [3] we consider the two-loop two-point functions, 
where we explicitly compute the wave-function renormalization and mass-shift. In Sec. H] 
we derive the near-flat limit of the conjectured S-matrix. In Sec. |5]we find the one-loop 
four-point amplitudes while in Sec. [6] we find these amplitudes at two loops. In these last 
two sections we also show that these results are in agreement with the results in Sec. HJ 
In Sec. [7] we present our conclusions. We also include several appendices which contain 
some of the technical details of our calculations. 



2 Near-flat-space model 

Our starting point is the relatively simple light-cone action for the reduced model of 
[14] (in the notation of [16]): 

+7 (Y 2 - Z 2 ) ((cLF) 2 + (d-Z) 2 ) + % 1 (Y 2 - Z 2 )tfjd.^ 

+i 7 ^(a_r i 'r' + d-tfr) (Y^r' - z i r i )ij) 

-^(ipr i,j 'ip ipr i,j 'ip - i)F ij i) ipr ij tp) . (2.1) 

The bosonic fields Z % and Y 1 ' correspond to transverse excitations in the AdS§ and S 5 
directions respectively and the fermions, ip, are Majorana-Weyl SO (8) spinors of positive 
chiralityQ The action in 12.11 is not invariant under world-sheet Lorentz transformations, 
but it is invariant under 8 independent linearly realized supersymmetries. 

This action is the same as the near-flat space truncation of [14], however as in [16], we 
have introduced the parameter 7 by rescaling the worldsheet coordinates and furthermore 
we have integrated out the half of the original sixteen fermions which occured only 
quadratically in the action. The near-flat space action action of [14] was obtained from 
AdSs x S 5 string sigma-model by expanding about a constant density solution boosted 



'"See App. [C]fbr a more complete description of the relevant conventions and notations. 
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with rapidity ~ A 1 / 4 in the a direction and so the above truncation should be equivalent 
to the full theory in the near-flat limit, 



provided we set 

and the mass, m, to be unity. 



p_~^A , P+~^j ( 2 - 2 ) 



" = 7t (2 ' 3) 



3 Two-loop propagator 

We now turn to the computation of the two-loop correction to the propagator. Firstly, we 
confirm that this leads to the expected mass shift and therefore the expected corrections 
to the dispersion relation. Secondly, for our two-loop scattering computation in Sec. El 
it is necessary to know the residue of the pole in the propagator, which we determine 
here as well. 

The dispersion relation in the original sigma model is expected to be 



mWl + — sin >\m z + p z -- — - . (3.1) 

7^ m V 3m 2 



The second expression is the predicted exact dispersion relation in the near-flat limit 
(12.21) . We will now derive this dispersion relation from a Feynman diagram computation 
in the model (12.11) . This computation shows for the first time the emergence of the sine 
in the dispersion relation from the perturbation expansion of the string sigma-model. 
The first correction to the propagator is of order r y 2 and the corresponding diagram 



is the sunset diagram drawn in Fig. 3(b) on page |2"U1 Doing the combinatorics for the 

(3.2) 
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bosonic and the fermionic propagator, respectively, leads to 

A b (p) = 32* 7 2 [2p 2 _ (2/ 110 + J 200 ) - p- (/ m + I 2W ) + (4/ 2 ii + /220 + / 
A f {p) = 16z 7 2 [p 2 _ /100 + 2p_/ 200 + 6/in + 14/2io + /300] 

where 

[ dPkdPq (k-Y (q-Y (p_ - k- - q-Y 

rst{P> J (2vr) 4 (k 2 -m 2 )(q 2 - m 2 )[(p - k - q) 2 - m 2 ] ' 1 ' ' 

This integral is the sunset diagram with r, s and t powers of the three momenta inserted 
into the numerator, cf. App. IB. 21 Some of these factors originate from derivative cou- 
plings, others are due to the extra power of p_ in the fermionic propagator. We can 
simplify the expression for the amplitudes using the identity 

P-Irst — Ir+l,s,t + Ir,s+l,t + Ir,s,t+1 ■ (3-4) 

Applying this identity repeatedly, we find that the amplitudes simplify to 

64 64 
A b (p) = y^7y-/ooo(p 2 ) , A f (p) = yHV-/ooo(p 2 ) , (3.5) 
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where Iqoo is a function of p 2 only. It is interesting to see how the very different structures 
in (13.21) reduce to essentially the same expression. We perform this integral in App. IB. 21 
and find for the on-shell amplitudes 

In order to find the corrected dispersion relation, we consider the iteration of sunset 
diagrams (13. 5ft . Via a geometric series this leads to the corrected propagator 

i i %Z(p ) 

p 2 - m 2 + fry 2 pi/ oo(p) 2p + - £(p_) 

where there is an extra factor of p_ in the numerator for the fermionic propagator. The 
right hand side of (13.71) defines the position 2?(p_) and the residue Z(p_) of the pole in 
the propagator in the 2p + plane. Note that for our definition of the light-cone momenta 
(1C.1I) . 2p + is the appropriate "energy" for time evolution in cr + direction. 

The dispersion relation is determined by the pole in the propagator. To order r y 2 we 
only need the on-shell value flB.7j) of the integral I 000 and find 



P+ (p_) = i £ ( P _) = ^-^. (3.8) 

Using e 2 — p 2 = 4p + p_, we convert this equation into the form e(p) and find that this 
dispersion relation exactly agrees with the prediction in (13. ip . 

For computing the residue we also need the on-shell value of the first derivative of 
iooo with respect to p 2 . Taking this integral from flB.8p . we find the wave-function 
renormalization to order r j 2 to be 



Zip-) 



2p_ 



7 2 



1 m 4 U 2 12 ' '' 



A 



(3.9) 



This correction is an important contribution to the two-loop amplitudes which we com- 
pute in Sec. El It will turn out that this correction cancels the entire wineglass contri- 
bution in the t-channel. 

4 SU(2|2) S-matrix 

The SU(2|2) scattering matrix is expressed in terms of the following kinematic variables^: 

1 + Vl + P 2 n . np 

x±(p) = b e ^ ' P = sm — . (4.1) 



2 We use the string normalization of momenta, which differs from the spin chain normalization in [5] 
by a factor of 2n/^/\. 
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For the S-matrix components, we use the conventions of [15]: 



■ 



Sg = A 6 C X + B 5 d a 5i , S« = C e^ 5 

SI 5 , = D 5% + E 5 S J} , fig = F ^e cd , (4.2) 

ScS s~i rc c5 ri'ld jj cc/r7 

a p = C6 a 6p, S a(j = Hd a d (3 , 

S£ = LS2S*, S£ = Ktf a SZ. 
The explicit expressions for matrix elements are 

A = 



B 



C 



D 



E 




i(x + — xJ)(x' + — x'_) 1 x' + — x + **pi 



T]T]'X-X_ 1 — — X_ — X + 

G = ^ + e s/x , H=-L^ -e , 

x'_ — x- rf x + - x- . A oX 

L = e vx , K — — , 4.3 

ar_ — x+ 77 x_ — x + 

where x± = x±(p), x'± = x±(p') and 

I I 1/2 / I / / 1 1/2 (a a\ 

7] = |x_ — x+| , 77 = |x_ — x + | . (4.4) 

The sigma-model scattering matrix is the tensor product of the two SU(2|2) S- 
matrices. The world-sheet scattering amplitudes are thus quadratic in the A, B, C,D, . . .. 



3 Comparison with the explicit tree- level calculations [15] shows that the scattering in the sigma-model 
is described by the SU(2|2) S-matrix in its canonical form [21] and should include phase factors e ±7rlpl / v ^ 
and e ± 7 "Pi/v / A j-j^ multiply the S-matrix elements in particular combinations. In other possible forms, 
which are related to the canonical S-matrix by state-dependent unitary transformations [21], e ±7ri Pi/V^ ; 
giTripj/x/X are replaced by arbitrary functions of p, p' [6] (for instance by 1 as in the original proposal [5]). 
It is interesting to note that in the near-flat-space limit the phase factors scale away and can be dropped 
altogether. 
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In addition the world-sheet scattering matrix contains an overall phase factor: 



1 - 



1 - 



x , — X- 



(4.5) 



where 8(p,p') is the dressing phase. For reader's convenience we have written the action 
of $ on all two-particle states in App. [D] in order to see which matrix elements govern 
which processes. 

The dressing phase has the following general form [19,22]: 



0(P,P') = rs x(xr,x' s ) 



(4.6) 



r,s=± 



The function x(x, y) is anti-symmetric in x and y and can be expanded in asymptotic 
power series in rr/^/X. We only need the first three orders of this expansion: 



xi x ,y) 



2tt 



x-y)[l- 
1 dt 



1 



In 1- 



nt 



In 



xy / 

1 - tfxy - (t 



xy 



;i - tfxy -(t + y) 



21 



'1 



7r xy + 1 



t) 2 xy- (t + x) 2 (1 -t) 2 xy 
x-y 



it - v) 1 



3V\ xy - 1 (x 2 - 1) (y 2 - 1) 



+ . . . 



(4.7) 



The first line is the AFS tree-level phase [19] , the second line is the HL one-loop correction 
[20] and the third line is taken from [7]. The integral in the one- loop phase can be 
expressed in terms of the dilogarithms, but for our purposes the integral representation 
is more convenient. The first and last lines are part of BHL's even phase, while the 
middle line makes up the entire odd phase. 

In the near-flat limit, the kinematic variables x± approach —1. However, the S-matrix 
contains many expressions of the form x r x' s — 1 or x r —x' s which vanish at x r = —1 = x' s . 
Plugging in —1 for x±, x± produces singularities and we need to keep the next term in 
the expansion: 

I 171 

x ± = -l ±-=p_ + ... . (4.8) 

P- VA 

The second and the third terms are small compared to one (they are of order (9(A~ 1//4 )) 
and should be omitted wherever —1 does not cancel. 
We thus get 



where the matrix elements are as in 

(H2) withS 



A = 


1 + 


m 


i 

-V- 


V- 
]/_ 


-p- 
+ P- 


D = 


1 - 


in 


y_ 


p'- 

p'_ 


-p- 
+ P- 


G = 


1 + 


in 








L = 


1 - 


in 

7t p - 


-V- 







B = _ E Ain P-P'- 



c 



y/X P' 2 ~ P- ' 
o • 3/2 13/2 

lin p_L pj 
s/X p'_ + P- 

o • 3/2 /3/2 

ti — K — —-= — , 

V A V- - P- 

(4.10) 

We should stress that these expressions are exact in the near-flat limit. For comparison 
to the two- loop calculation in Sec. El we need to further expand in npi/\/X. 

When expanding the phase in n/^/X it is important to remember that it implicitly 
depends on A through x±, apart from the explicit dependence manifest in (14.71) . In 
particular the tree-level term in (14.71) contains a two-loop correction to the phase. The 
substitution of (14. 8 ft into (14. 6L (14.71) yields after a lengthy but straightforward calculation: 

af a 2tt , pL-P- 4tt 3 3 /3 (p'_ ~P-)(p'i +p-p'_ +p 2 _) 
V{P,P) = -^fP-P--j— t ^T^P-P 



,/X p'_ + p- 3A 3 / 2 (pL + p_ ) 3 

87T p 3 p' 3 ( p'i + pi p'_ 

+ T~ ~ 2 1 _ 12 2 ln 

\ p± — p± \ p'± — pi P- 

^ n 3 r/ 3 (P- ~ P-W- + ^P-P'- + P-) u x 

+ 3A 3/2^- (p'_+p_)3 l4 - iiJ 

Omitting (9 ^(-^=) 6 j terms this can be written in the following nice form, suggested by 
the main scattering term, 

, f n 1 + \fxP-P'~ P trt , svr pV 3 / p' 2 +p 2 p' LA 

0(p p =-zlll ; ^^ + ~r^2 2~ ( /2 r 111 — ' 412 

1 _ j* p v > p -- p - A p'i - pi V P /2 -pi P-J 
Vx p-+p- 

where the first (second) term comes from the even (odd) phased Equation (14.91) becomes 

3 /3 / 12 . 2 
8vr» P-P - I 1 P -+*■- j i 

a p'^-p-i V p'^ -p 2 ^ ; 

$ = S S®S with S = — 5— ( 4 - 13 ) 



t 2 ^2^/2 Z' 



2 

1 + ^pip'i 

where S is given by ( 14. 2ft . ( 14. 10D . At the end, the dressing phase almost completely cancels 
the main scattering phase, and the two-loop prediction for the scattering amplitude turns 



4 We chose to pull out a common factor of N — -^=p_p'_ p ~ ^ j from S. 

5 We suspect that the even part of the phase in 14.121 is valid to all orders in 7. We checked this 
by taking the near-flat limit of the BHL phase to order 7 . Furthermore, one can readily see that it 
solves the near flat limit of the even crossing relation of (2.13) in [7] . It should be possible to prove (or 
disprove) this fact by inspecting the integral representation of the phase found in [23] . 
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out to be rather compact. We should stress that (14.131) is only accurate up to 0{\/\ 2 ) 
(the full expression is expected to contain dilogarithms from the odd phase) while matrix 
elements (I4.10p are exact in the near-flat sigma-model. 

In order to facilitate the comparison with the results from the world-sheet compu- 
tation, let us discuss the first few orders of (14.131) . The ra-th loop contribution to the 
two-particle S-matrix is of order 7 n+1 = (-^=) n+1 and we denote it by §( n \ Now, we 
observe that the prefactor So in (" 14. 131) does not have a term of order 7 = and that 
the coefficients in (14.101) stop at order 7 = Hence, the tree-level contribution to the 

S-matrix originates only from the matrix elements in (I4.10p . the one-loop contribu- 
tion S^ 1 ) receives additional terms from the prefactor So and the two-loop contribution 
is of the form 



S (2) 



7T 
T 



2 12 

-pt.pl 



p_ + p- 
p'_ — p- 



8i p 3 _p'3 

IT p't. — 



p/2 _|_ ^2 ^ 

m — 



p/2 _ p2_ 



P- 



S (0) , (4.14) 



i.e. the two-loop piece reproduces the tree-level S-matrix multiplied by a factor that is 
universal for all scattering processes. 

We close this section by noting that the S-matrix can be put into a form that looks 
almost relativistic. Under boosts the momenta, derivatives and fields transform as 



P± 



a Tl p± 



ip — > yaip 



(4.15) 



where a is the boost parameter. If these transformation are accompanied by a rescaling 
of the coupling 7 — > a" 2 7, then the Lagrangian (12.11) is invariant under these transfor- 
mations. As a consequence the S-matrix can be written as a function of a momentum 
dependent, but boost invariant coupling 

7 = l(p,p') = 1P-P- 



and the relative rapidity 9 = arcsinh 
we find 



with 



S 



JL 

rn 



arcsinh ■ 



e t 



sinh ( 



1 + 7 2 coth 2 f 



(4.16) 

In^k Rewriting (T4~T3l and fOHjl . 
S®S, (4.17) 



A = 1 + ^tanh | , 
D — 1 — Z7tanh | , 
G = 1 + ^7 , 
L — 1 — ry . 



B = -E = 2^7 csch § , 
C = F = ry sech | , 
H = K = ij csch I , 



(4.18) 



It would be interesting to see if a proof of crossing symmetry can be obtained given this 
relatively mild breaking of the Lorentz invariance. 
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5 One-loop amplitudes 



In this section we present the general bosonic one-loop amplitudes and S-matrices for 
2^2 magnon scattering in the near-flat limit. These results generalize the case of 
ZY — > ZY presented in [16]. 

For all processes, there are three basic diagrams which are shown in Fig. [U We call 
these graphs the s, t and w-channel graphs. Within each of these graphs, there can be 
several contributions to the complete loop in that channel. However, summing over the 
contributions will lead to three basic structures for the one-loop amplitudes. The first 
of these is a structure associated with forward scattering, the second is a permutation 
structure and the third is a trace like structure. The latter two structures are related to 
each other through crossing symmetry. 

The one- loop amplitudes are relatively straightforward to carry out. For an amplitude 
of forward scattering type (for example Z\(p)Y\{p') — > Zx(p)Yi(p')) , the amplitude is 
found to be 



A 



(i) 

forward 



iP-,P'-) 



(p'_ (p'_ +p_)I 00 (p, -p') - 8p_p_I 00 (p,p) 



+ (P- - P-Y(P'- + P-)Ioo(p,p') 



(5.1) 



where Ioo{p,p') is the u-channel loop integral defined in (IB. II) . The s and t channel inte- 
grals are given by analytically continuing p'_ to — p'_ and letting p'_ — > p_, respectively. 
Substituting the results for the integrals into (15.1 1) gives 



-^■forward (P~ ' P- 



lQiY 



- 47 



pi-P'l 



P 2 -) 



p' 2 + pi , pL 

1 ^ ft In — 



p/2 _ ^2 



7T{P- 

2 P-p'_ (p'_ + P- ) (p 2 _ + p'l 
p'_ — p_ 



P- 



(5.2) 



This result was previously derived in [16]. 

The next type of scattering process is of the permutation type, where the outgoing p 
and p' are exchanged with a forward scattering process. In this case, summing over the 
contributions to the Feynman diagrams, we find 

•A§Lx(P-M = -!67 2 [p-P'-(P- + P-) 2 Ioo(P, -P') + P-P_(P- + P-) 2 ho(p,p')] , (5.3) 

where in these processes the contribution to the t-channel cancels out and the -u-channel 
integral comes with the same kinematic factor as the s-channel. Substituting for the 
integrals into (15.31) we arrive at 



•^pcrm (P- 



>P- 



*7 



2 P-P' 2 (p'_ + P- 



p_ — p^ 



(5.4) 



Finally the processes of trace type, which are of the form AA — ► BB, where A and 
B are any one of the fields and A and B are there conjugates is given by 



AlLe{P-,P~) = !67 2 P-P-(P- -P-) 2 ho{p,-p) +p-p'^{pi -P-)' 2 Ioo(p,p') 



(5.5) 



10 



p P P P P P 




p' p' p' p' p' p' 

(a) s-channel (b) t-channel (c) u-channel 



Figure 1: One-loop graphs. The basic structure for s, t and u channel graphs. 



which after substituting for the integrals gives 

X'L( P - P -)=8 7 ^- P -, ( ^- p -> . (5.6) 

The amplitudes in (15.41) and (15. are related by crossing symmetry by taking p'_ — > 
—p'_. However, there is a subtlety in the analytic continuation, since the amplitudes 
were obtained by continuing around a log cut. When continuing, say, A[]l ce (p-,p'_) to 

Afrlce(p-, —p'-), one continues onto a different branch, hence leading to an extra minus 
sign. 

6 Two-loop amplitudes 

In this section we compute the two-loop amplitudes for various four-point processes 
and show that there is complete agreement with the S-matrix results in Sec. HI One 
consequence of the structure of the S-matrix is that the two-loop amplitudes should be 
related to the tree amplitudes by a universal factor / j 2 F(p_,p'_), which we will explicitly 
show. 

In order to obtain the S-matrix, one must take into account the wave-function renor- 
malization of the external legs as well as a Jacobian factor that arises when converting 
(^-functions for overall conservation of energy and momentum to ^-functions for indi- 
vidual momenta. Moreover, there is a two-loop contribution to this Jacobian due to 
the two-loop mass-shift. The contributions from the wave- function renormalization and 
Jacobian will cancel off against certain terms in the amplitude to give very compact 
expressions for the S-matrix. 

Since all interaction terms in (12. ip are four-point, the general structure for the two- 
loop Feynman diagrams have the form shown in Fig. [2J The diagrams fall into the 3 
general classes, "double bubble" , "wineglass" and "inverse wineglass" for each of the s, 
t and u channels. The bosonic vertices all come with two powers of p_, a vertex with 
two bosons and two fermions has one power of p- , while the four-fermion vertex has no 
powers of momenta. The fermion propagator also comes with a factor of p_, therefore 
the amplitudes will all have world-sheet spin —6. Naive power counting might indicate 
that these diagrams are divergent, however the two-dimensional Lorentz invariance of 
the free theory insures that these divergences are not there. 
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(a) (b) (c) 




(d) (e) (f) 




Figure 2: Two-loop graphs. The first line shows the s-channel diagrams, the second line the 
t-channel and the last u-channel. To the diagrams in the first column, we refer to as "double 
bubble" , to the second column as "wineglass" and to the third as "inverse wineglass" . 
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Let us start with the easiest set of diagrams to evaluate, the t-channel double bub- 
ble. For these diagrams, no external momentum flows through the internal propagators. 
One can argue that there must be at least two powers of the internal momenta in the 
numerators of the two-loop integrals, which by the Lorentz invariance of the free theory, 
must be zero, and so the t-channel bubbles all have -4. bubble (p,p') = 0. 

The next set of diagrams we consider are the u and s channel double bubbles. Differ- 
ent processes have different combinatoric factors contributing to the loop integrals, but 
their final results all reduce to the same form, with the w-channel given by 

A h u uhblc {p,p') = j 2 A^(p-,p'_) F bubhlc {p_,p'_) , (6.1) 

where A^°\p-,p'_) is the tree-level amplitude for the corresponding process and 



^bubble 



(P~,P- 



[2(p_ -p^) 2 I 00 (p^,p'_ 



(6.2) 



I rs (p-,p'_) is the one-loop w-channel integral defined in ( IB. II) . The s-channel can be 
obtained easily from the w-channel result by continuing p'_ — > —p'_ in F (p_,p'_), 
but not in A^°'(p-,p'__), resulting in 



A 



bubble 



(p_,pL) = 1 2 A^{p-,p'J)F h ^{p^-p>_) . 



(6.3) 



Combining the double bubbles together and substituting the expression for I 0Q (p_,p'_ 
in (IB.2p . results in 



A 



bubbles 
stu 



7 2 i (0) (p-,/ 



p_p' 2 (p'_ + P- ) 2 2ip 2 _p' 2 (p'_ 



V- 



(p'_ -p-) 2 



ic(p'_ — p_) 2 



P- 



+ 



2plp' 2 (pi + 6p 2 _ p' 2 + J9 74 ) , 2 p'_ 

— In — 
V- 



7l 2 (p 12 — p 2 _) 2 



(6.4) 



The wineglass diagrams are computationally more challenging because their loop 
integrals do not factorize into products of one-loop integrals. Nevertheless we are able 
to obtain compact expressions for these as well. Like the double bubble diagrams, all 
2 — > 2 processes have the same proportionality factor to their tree level amplitude. For 
the ^-channel wineglass, we find the expression 



Ar e (p, P ') = l^ 2 A^\p-MF^{p^p'_ 



(6.5) 



where 



F wine (p-,p_) = 16 



^p^p' 2 Wo (p- ,p'_) + 8p-p'_ (j>'_ + p^ ) W\ (p_ , p'_ ) 
(p 2 _ + 6p_p'_ + p' 2 ) W 2 (p-,p'_) 



(6.6) 



The wineglass integrals W r (p-,p'_) are defined and discussed in App. IB. 31 Different 
processes have very different combinations to reach this same final form in (16.51) and 
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( 16. 6ft . The s-channel wineglass is again related to the n-channel form by analytically 
continuing p' — > —p' in F wme (p_,p'_), 



AT ne (p,p') = -fA^fr-M F win >_, -p'_ 



(6.7) 



Likewise, we also find that the t-channel wineglass has a simple relation to the other 
wineglass diagrams, namely we simply set p'_ = p^ in F wmc (p_,p'_), giving us 



ATM) =7 2 A {0 \p^p'„)F™ c (p„,p_) . 



(G.i 



For the inverse wineglass diagrams, it is straightforward to show by the symmetries 
in the diagrams that 

A i r Tsc (p-, P '-) = Ar c (p-,p'-), 

Ar ise (P-,P-) = AT lc (p-,pi) , (6.9) 

while the t-channel inverse wineglass is 

Ar CTSC (p,p) = j 2 A(°\p-,p'„) F™ c (p'_,p'_) . (6.10) 

Putting together the terms in (16. 5p . (16.71) and (16.91) and also using (16.61) and the 
expressions for W r {p-,p'_) in App. IB. 31 , we obtain the combined su wineglass 



A 



wineglasses 
su 



7 2 AW(p„,p'_) 



2 12 
PJPl 



ip z _p'^ 2ip 2 _p' 2 (pt + Qp_p + p ) p_ 

H — — m — 

P- 



7r{p' 2 — p 2 ) 7T(p' 2 — p 2 ) 2 

/4 



2p 2 _p' 2 (pt + 6p 2 _p' 2 + p'Z) 2 V 

2^2 _ n 2 \2 111 ~ 



ir 2 (p' 2 — p 2 _) 



Combining the t-channel wineglass with its inverse gives 



^winegla 



7 2 ^(p-,^; 



7T- 



1 

12 



(p'*+p 4 



(6.11) 



(6.12) 



and then combining this with ( 16.41) and ( 16. lip , we reach the final two-loop amplitude 

A^{p^ P '_) = ^ bics (p_, P '_) +A:r slasses (p-,p'-) +A7^ iasscs (p^p'_) 

p 2 _pV(p'_ +p_) 2 



7 2 ^°)(p_,p'_) 



(p'_ — p_ 



+ 



■tip^p' 3 



7T(p' 2 — p 2 _) \ p' 2 — p 



p' 2 + p 2 - , p'_ 

1 ^ — m — 

P- 



2 .'2 I * \ i 'i I 

-P-P- + I ~ yjj ) (P- + P- 



(6.13) 



One should immediately note that the (In — ) 2 terms that appear in (16 .4p and (16. lip , 
but which are absent in the two-loop S-matrix in (14. 10114. 13p have canceled off in the final 
amplitude! One can also easily see that the first line of (I6.13P has precisely the right 
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form as ( 14.14ft . The first term in the second line is accounted for by a Jacobian factor, 
while the second term in this line, which is due entirely to the t-channel contributions, 
is compensated by wave-function renormalization of the external legs. In fact, the renor- 
malization of the legs with momentum p through them cancels off with the t-wineglass, 
while the renormalization of the p 1 legs cancels against the inverse t-wineglass. 

The Jacobian arises because the amplitudes come with factors of 6 2 (P^ nt — P£), while 
the S-matrix is written with factors of <5(p_ — qj)5{p'_ — q'J). These are related by 



1 / dp' + dp-\ 



out ± out/ 



S(p_ -q_)5(p'_ - q'J) . 



(6.14) 



2 \dp'_ dp 

Taking into account the two- loop dispersion relation in (I3.8p . we find for the Jacobian 



1 / dp' + dp^ 



2 \dp'_ dp. 
The full S-matrix has the form 



2 m l2 



2p±p 



m 2 (j/2 _ ^2 



1,7 2/2 

1 + —;P-P- 
m 4 



(6.15) 



S = 1 



1 / dp', dp^ 



2 \dp'_ dp- 



-i 



Z(p.)Z(p'_)A. 



(6.16) 



Thus, after setting m = 1 and substituting in (I6.15P and (I3.9p . the two- loop contribution 
to the S-matrix is 



S (2) 



P-P_ 



2(p' 2 -p 2 _ 

Using the result for in (I6.13[) . we reach the final expression 

S (2) =7 2^(0 ):r _^ 



(6.17) 



•^3 

2(p' 2 -p 2 _) 



p_ — p- 



P-+P- \ ^ P-P 



7T p' 2 — p 2 _ 



p'1 + pL , p_ 

1 ^ 7T m — 



p' 2 — p 2 - 



V- 



(6.18) 



which agrees precisely with the conjectured form (14.141) . since 

Sj(o) =^(o)_^ 



V-P- 



7 Conclusions and outlook 

The sigma-model describing the super-string on AdSs x S 5 is a rather complicated theory 
and calculating the complete quantum S-matrix remains a formidable problem. Fortu- 
nately consideration of the near-flat space limit, as described in [14], results in significant 
simplifications which make loop calculations feasible. The reduced sigma-model has at 
most quartic interactions and the right movers essentially decouple from the interacting 
left-movers. Just as for the full string theory in the light-cone gauge the reduced model 
is not Lorentz invariant however if boosts are combined with a rescaling of the loop 
parameter the action is indeed invariant. This can be seen in the world-sheet S-matrix 
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which depends only on the difference of rapidities and an effective, momentum depen- 
dent, coupling. Furthermore the simplified theory possesses at least (0,2) worldsheet 
super symmetry. 

As an important step in the calculation of the S-matrix we computed the two-loop 
two point function with the corresponding mass shift and wavefunction renormalisation. 
This is an interesting result in its own right as we can explicitly see the modification of the 
relativistic dispersion by the sine function at higher powers of the momenta. In the gauge 
theory description the sine function arises naturally from the intrinsic discreteness of the 
spin-chain and indeed from the point of view of soliton description [18] the momentum is a 
periodic variable as it corresponds to the angular separation of the string endpoints. This 
is however the first case where the sine function has been seen to originate from quantum 
corrections to excitations about a plane-wave vacuum. Additionally in calculating the 
full SU(2|2) S-matrix we are able to check that the symmetries of the classical theory 
are realized at higher loop order. 

Given the central role of the world-sheet S-matrix in recent developments of our 
understanding of the AdS/CFT correspondence it certainly interesting to extract as 
much information and intuition from this reduced model as possible. The spectacular 
agreement of our calculations with the appropriate limit of the conjectured exact S- 
matrix of [7], [8] provides further strong evidence in favor of its validity. It should be 
straightforward, though perhaps technically challenging, to extend the loop calculation 
to even higher orders which would provide yet further confirmation of the complete 
S-matrix. However, given that the theory is presumably integrable, it may be more 
profitable to try to find a complete solution using more non-perturbative techniques 
perhaps along the lines discussed in [24] . This would allow one to answer an outstanding 
issue not addressed by the perturbative calculation, that of the pole structure of the 
S-matrix. Although we consider the near-flat space limit which interpolates between 
the plane wave limit and the giant magnon regime we do not see the double poles of 
the S-matrix corresponding to exchange of BPS magnons [23]; which would require a 
resummation of the entire perturbative expansion. 

Note added While this paper was being prepared for publication we received [25] 
where the study of two-loop quantum corrections to the energies of classical string solu- 
tions was initiated. 
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A S-matrix elements 



A.l Bosons 

We write the action of the T-matrix, which is defined as S = 1 + T, onto all bosonic 
initial states. We omit fermions in the final states. Using an SO (4) ® SO (4) notation, 
we define the matrix elements as follows: 

T\Z,Z' 3 ) = \ZiZ$I zz + \Z 3 Z[)P ZZ + 5 i3 \Z k Z' k )T zz + h 3 \Y k ,Y' k )K zz , (A.l) 

T\Z,y d ) = \ZiY;,)I ZY + \Y f Z[)P ZY , (A.2) 

T\Y V Z' 3 ) = \Yi>Z' 3 )I YZ + \Z 3 Yl,)P YZ , (A.3) 

T\Y,Y^) = \Y V Y;,)I YY + \Y 3 ,Y[)P YY + b V3 ,\Y k Yl)T YY + 5 Vjl \Z k Z' k )K YY (A.4) 

The world-sheet computation yields 
r(o) _ _ 7 (o) _ _ 2 . P-Pl(pg + P 2 -) 

A ZZ — YY — Z H p ,2 _ p 2 

p (o) _ _ p (o) _ _ 4 P-P' 2 

r ZZ — YY — 12 _ p 2 

Tzz = = +^7 a _ ^p, 

Kfz = -KPy = 



(i) _ 9 2 P-P' 2 (p' 2 +P-) 8^7 2 pip' 3 f _ p a + p 2 p'_ 

zz ~~ "r-'yy ~~ ^7 / / \ 9 + /2 2 I /2 2 

[P_ — P-Y IX p_—p Z _ \ pl—pt P- 



p(D _ , p(l) _ P-P- 

^ 1 2 = +4 1 , = +47 2 



2 



4 1 i = +4 1 ]- = -47 2 ^ ? 



(p'_ + P-) 2 
(p'_ + P-) 2 



(2) r(2 ) , . s P^g^.+p^^+p 2 ) ^pV-V-g + p 2 ) A p'j + p 2 , p 

ZZ- *YY-+*1 {p'_-p_Y + 7T (p'_ 2 -p 2 ) 2 V p' 2 - P 2 _ n p 



,(2)_ p (2) _ ^- sP-P-(p'-+P-) , 32 7 3 pV- p'_ 2 +p 2 ^p'_ 



Pzz " P nr - +4*7 " (p 7 _" p _ )3 + ^ ^ " p2 )2 ^1 ^ _ ^ 2 In ^ 

T (2) _ _ T (2) ,. 3 pip-(p'- + P-) 32 7 3 pVj / p'j+p 2 p^ 

ZZ " 7 (p'_"P-) 3 7T (p'_ 2 -pl) 2 V P- 

Kgl = -K& = 
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4°r = -41 = -2ijp-p'„ 

p(0) _ _ p (o) _ n 
r ZY — r zy — u 



r(l) ,r(l) o 2 P-P'-jP- + Pi) , 8^7 2 / P^ + pj , pi 

izy ~ +iyz ~ ^ 7 - P _y + tt p'„ 2 - pi v p- 2 - p 2 - P- 

p(D _ ,p(D _ _ 4 ^2 P-P' 3 

r (2) _ r (2) , 9 - 3 P-P'-(P- +P~Y , 167 3 plp_ A P^+p 2 , P^_ 

Zy ~ YZ ~ + 11 ( P '_ - J9_) 2 + 71 p'l - pi V P' 2 ~ Pi P- 



p(2) _ _ p (2) _ n 

-'zy — r yz — u 



These coefficients have to be compared to the S-matrix elements (14. 13ft in the follow way: 

(A.5) 



1 + Iyy 


= S A(A + B) , 


l + Izz 


= S D(D + E) 


Pyy 


= S B(A + B) , 


Pzz 


= S E(D + E) 


T Y y 


= S AB, 


Tzz 


= S DE, 


Kyy 


= —Sq C 2 , 


Kzz 


= —S F 2 , 


I + IZY 


= S L 2 , 


1 + Iyz 


= SqG 2 , 


PzY 


= S K 2 , 


Pyz 


= SqH 2 . 



Here So denotes the pref actor in (14.91) . We find perfect agreement. Note that we are 
sensitive to all matrix elements, even though we concentrate onto the scattering among 
bosons. This is because the field Z actually carries two fermionic indices in the su(2|2) 2 
notation. 



A.2 SU(2|2) subsector 

We now extend our considerations to include processes involving fermions, however, we 
restrict ourselves to a single su(2|2) sector. Granting the group factorization of the full 
S-matrix, this is a sufficient test of the supersymmetries at higher loop orders. 

As described in App. [Cj we identify in the worldsheet theory the fields 4> a and Xa 
spanning an su(2|2) sector. We calculate the matrix elements defined as follows: 

$|0 a 0' 6 ) = SoAM) + SoB\<p b <p' a ) + SoCSabS^lXaXp) , (A.6) 

$\XaXp) = SoV\ X aXp) + So£\xpX a ) + So^e^e *^) , (A.7) 

SI&XJ) = SoGl^aXp) + SoH\xM , (A.8) 

S|x«$) = SofC\<j)bX a ) + S C\x a (j>' b ) ■ (A.9) 

For the sake of brevity we will not record the individual amplitudes but simply state 
the final results for the .S-matrix elements, noting that they agree with the all-order 
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prediction from the dual spin chain description. There is a common contribution to each 
element of the form So = 1 + iS where 

71 p_ — pt \ P_ — pt P- J 

in addition to the individual contributions 

P-p'_ {p 2 _ + p' 2 ) n p 2 _p' 2 {p 2 _ + p'l) o P-pf-ivt + P-P 1 - + P-p'l + p' 3 ) 

A = 1 + 2«7 -= = 2 7 — — — — + 2tY 



p'l — p 2 _ (p'- — P-) 2 (p- _ p~y 

„ V 2 V' 2 n P 3 P' 3 n P A P' A {P- + V' ) 

B = -E = 4 ?7 - 4 7 2 - 4z 7 3 ^ 

pz—pz. \P-~P-) \P- ~ P-) 



C = ^ = 2^7 -r-j 2 7 7^ r - 2^7 j-, t~ 2 

p_+p_ [P- ~ P-) \P- — P-) 

^ , . A . P-P- , 2 P-P- A . 3 P-P-(P- + P-) 

V = 1 + 4^7 — T - 4 7 — - 4?7 — — - 

pi — pi [p- — P-) \P- — P-) 

r , . n . P-P- _ 2 P-P-(P-+P'-) o- 3P-P-(P- +P-Y 

g = I + 227- 2 7 — — ; 2«7 — — — - 

p'_ _ p_ _ p_ y _ 

n, r n . iP-P-)^ . 2 - 3 (P-P-)Hp- + p'-Y 
U = K. = 2?7 + 2 7 — 2iY 7-, t- 3 

p_ — P- \P-~P-) \P- — P-) 

r , . n . P-P- 2 P-P-(P- + P-) n . 3 P-P-(P- +P-f 

C = 1 + 2^7 2 7 — — ; 2«7 — — — - . 

p_ — p~ [p_ — P-) \P- — P-) 

In considering a single SU(2|2) sector the full S-matrix is 

$ = S S® (A + B) (A.ll) 

as one index in the tensor product is kept fixed by the scattering. Thus we can write 
these elements in a simple compact form in terms of the S-matrix defined in 14.101 

$ = rr S g> 1 , (A.12) 

1 -l-1P-P_-T 



P--P- 



and which of course is in agreement with the AdS/CFT prediction to this order. Thus we 
see that the symmetries are preserved to at least two-loops in the reduced sigma model. 



B Integrals 



B.l Bubble integral 



We consider the bubble integral, cf. Fig. 3(a), for two inflowing momenta p and —p' as 
appropriate for ^-channel processes. With r and s powers of momentum inserted, the 
integral reads 

d 2 k 



(k-) r (p_ — p'_ — k-Y 



(27r) 2 (k 2 -m 2 )[(p-p'-k) s 



m 



(b.i; 
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k, r 




(a) bubble (b) sunset (c) wineglass 



Figure 3: Loop diagrams 



These momenta originate from derivative couplings and fermionic propagators. However, 
it turns out that all amplitudes simplify such that we only need to explicitely compute 
ioo> which is immediately found to be 



t ( '\ i P-P- i P- tuo\ 

ho{P,p) = 7 ) 2~2 2" ln— • ( B - 2 ) 

lum 1 p_ — p_ p_ 

In the s-channel, the inflowing momenta are p and p'. The integral 1qo(p, —p') is obtained 
from ( IB. 21) by analytically continuing the logarithm. In the t-channel, the total inflowing 
momentum is zero and we obtain from ( IB. 21) in the limit p' — > p: 

^p) = sss? • < B - 3 > 



B.2 Sunset integral 



The general sunset diagram, Fig. 3(b), is defined as 



rst[P) J (2vr) 4 (k 2 -m 2 )(q 2 -m 2 )((p-k-q) 2 -m 2 ) ' K ' } 
There is the relation 

P-Irst = Ir+l,s,t + Ir,s+l,t + Ir,s,t+1 (B-5) 

between different integrals which follows immediately from taking the p- on the left hand 
side into the integrand and writing it as + g_ + (p_ — A;_ — Using this identity, it 
is possible to reduce all sums of sunset diagrams that occur in the two-loop propagator 
to Jqoo- We solve this integral by introducing three Feynman parameters 

, If 1 dxi dxo dx 3 S(xi + x 2 + x-x — 1) ,„ , 

WP 2 = T7— 2 / 2 ( \ , \ 1 — • B.6 

Observe that this integral depends only on p 2 . On-shell the value of the integral is 

WP 2 = - 2 ) = ^ • (B.7) 
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Apart form this, we also need the on-shell value of the first derivative of Jooo with repect 
to its argument, which is given by 



j ooo(p 2 = ™ 2 ) 



1 dx\ dx 2 dx 3 8(xi + x 2 + x 3 — 1) Xix 2 x 3 



(47r) 2 m 4 J [xix 2 + X1X3 + x 2 x 3 - x 1 x 2 x 3 ) 2 
3/1 1 



64m 4 V 7T 2 12 



(B.8) 



B.3 Wineglass integral 



The wineglass diagram as drawn in Fig. 3(c) reads 

d 2 kd 2 q (k^) ni (q^) n2 (k„ + g_ -p^) n3 (p'_ - k- -q-] 



W rstu (p,p') 



(2tt) 4 (k 2 - m 2 ) (q 2 - m 2 ) [(k + q - p) 2 - m 2 ] [(p' - k - q) 2 - m 2 ] ' 

(B.9) 

We note the identities 

W rstu (p,p') = W srtu (p,p') = (-l) t+u W rsut (p',p) = (-l) r+s+t+u W rstu (-p, -p') . (B.10) 

All sums of wineglass integrals that occur in the two-loop amplitudes can be reduced 
to combinations of the following three terms which we compute by standard means and 
find 



W (p,p') = Woooo(p,p') 



p^p_ 
167r 2 m 4 

W 1 (p,p')=W 1000 (p,p' 
P-p'_ 



7T 



+ 



In* 



P^P_ 



_4(p'_+p_) 2 {p' 2 —p 2 ) P- (p' 2 — P-) 2 P-_ 



167r 2 m 4 



7T' 



P-P^ 



111" 



P. 



8(p'_ + p^) 2(p'_ — p){p' 2 — p 2 _ ) p 



W 2 {p,p') = W uoo (p,p') + W 20Q0 {p,p') 



P-P_ 
167r 2 m 4 



7T 



p^p_ 



12 2(p'_-p-] 



P-. 



(B.ll) 



(B.12) 



(B.13) 



C Notations 

In this section we summarize several of the notations used throughout the main text and 
record several useful results. We make use of the light-cone coordinates and momenta 

a ± = a°±a 1 , p ± = ±(p ± Pl ) (C.l) 

so that the worldsheet metric is rj^ = (+,—). We also use the notation e = po and 
p = pi, and bold-face for world-sheet two- vectors like p = (po,Pi)- 
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It is convenient to perform quantization in world-sheet light-cone coordinates with 
o + as time and where the target space fields have the mode expansions 

dp- 1 



Ziitr) 
Y v {tr) 



2tt 

dp- 1 

27T y^p~ 

dp^ 1 



o»(p-) e" 



rper 



Oi(P-)e" 



a i ,(p_)e- ip " T + 4(p. 



ip<T 



tp<T 



ip<T 



2tt 

The free bosonic and fermionic propagators are 

i ip 



P- 



-ipcr 



p2 _ m 2 



(C.2) 
(C.3) 
(C.4) 

(C.5) 



and the free dispersion relation is 2p H 



2p-" 



We use the following representation for the 16 x 16 7-matrices 



f 1 = e X e X e x e 
f 2 = 1 x Ti x e x e 
r 3 = lxr 3 xexe 
T 4 = ri x e x 1 x e 



r 5 = r 3 xexlxe 
r 6 = exlxrixe 
F 7 = f x 1 x t 3 x e 

T 8 = 1 X 1 X 1 X T! 



(C.6) 



with 



1 
-1 



1 

1 



T 3 



1 
-1 



(C.7) 



We also define r 9 = J 11 / 12 ■ ■ ■ T 8 and Pl,r = |(1 ±-T 9 ). The fermion ^ is a real, positive 
chirality spinor and hence has eight real degrees of freedom. 

The su(2|2) sector considered in Sec. IA.2I is spanned by the bosonic fields 



T2^ + m 



(at 



and the fermionic fields, Xa, which are most easily defined in terms of the projection 
operators, 

1 

2 
1 

2 

1 



PL 



in 



(i±r x r 2 r 3 r A ) 



1 ± - (r 12 + r 34 ) 



(C.9) 



such that, 



Xs 



piiipIipi_p L ^ . 



(CIO) 



These correspond to the fields (fi a = Y a i and Xa = m the notation of [15]. 
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D S-matrix action 



We spell out the action of the S-matrix onto the entire set of two-particle states in su(2) 4 
notation, cf. [16]. This serves as a reference for which processes can occur. Taking into 
account that the coefficients B, E, C, F, H, K are of order 7 = ^S=, we see that some 
of the processes are absent at tree-level. The terms that are present at tree-level are 
printed in bold face. To simplify the formulas, we suppress the So that multiplies all 
right hand sides in the following. 

Boson-Boson 

$1^) = + # \Y ah Yl h ) + AB \Y bh Y^) + AB \Y ab YQ + B 2 \Y bb Y^) 
+ AC e ab e^\Y^T' sb ) + BCe ab e^\T^) 

+ ^ £ /|^> + BCe^\^'J 
+ C 2 e ab e, h e^\Z^Z > 5S ) 

S\Z aa Z'^} = + D 2 \Z aA Z^) + DE\Z^Z'J + DE\Z a$ Z^) + E 2 \Z p$ Z' a „) 
+ DFe a ^ d \Y a ,Y^) + EFe^ d \T Pc T'J 
+ DFe^l*^*^) + EFe a pe cd \V c$ V' da ) 

+ F 2 e aP e^ d \Y ct Y' dd ) 

$\Y aa Z>^) = + G 2 \Y ah Z'^) + GH \TpP'J + GH \^ a0 Y^.) + H 2 \Z^Y' ah ) 

S|Z Qd ^> = +L 2 \ Z a< yY'-) + KL |flWTM + KL \Y ab ^) + K 2 \Y bb Z' adt ) 
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Fermion-Fermion 

S|SW?y = +AD \# a& #^) + BD \# b „*'J + |fl^ff^> + BE |^ d > 

+ AFc^iy^) + BF e cd ^\Y bc Y'J 
+ CDe^e^Z^Z^) + G£e ab e^|Z 7/ ^ d ) 

+ GFe a& e^e^|r 7< ^) 

s|r aa r; b ) = + ad \r ah r^ b ) + al |r^r^.) + bd \r ab r^) + 5£ \r pb v aa ) 

+ CD e ab ^ s | Z^ZV > + CE e hb e^ \ Z^Z'J 
+ CFe hb e a(3 e cd e^\^'J 

s|^ a r; b ) = + gl |^r;.) + el \z^y>j + gk \Y ab z>^) + m |r^ d > 

S|r Qd ^> = + GL |T^.) + GK \Y bh Z'J + iJL \Z^YQ + M |S^7^> 

Boson- Fermion 

S\Y aa ^} = + AG |F a ^) + £G \Y b& *' a$ ) + AiJ |^F4> + BH 
+ CGe^r^Z'^) + Gtf e^lZ^) 

$|^> = + AG |F a ^.) + AH \Yp h Y^) + BG > + BH \T^.) 

+ CGe^'^Z'^ + CHe. b e^\Z^'J 

%\*a^ b ) = + AL \# adi Y> b ) + fiL + AK \Y ab *> b& ) + \Y bb V>.) 

+ GL e ab e^ | Z^T' 6h ) + GK e a6 e^ jT^Z^) 

S|T Qd ^) = + AL \T ah Yl b ) + AK \Y bh T'J + \T ab Y> h ) + 
+ CLe. b e^\Z a ^') + GKe^^Z^) 
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S|iWO = +DL \Z a6 &>.) + M \# b «Z'.) + EL\Z a ^' b& ) + EK \^Z' aa ) 



s|z aQ ri fe > = + |z ad r'.) + fl |z^r' .) + dk \r ah z') + f^ \r 0b z> aa ) 



S^Zp = + -DC |#a*Z^> + DH\Zp^'J + FG |a^Z^> + FF |Z^< d > 
+ FGe.^\Y a6 r^) + FHe^ d \T, c Y'J 

%\r a& z' 0$ ) = +dg \r ah z'^) + eg \r^z'j + dh \z a0 r^ > + ff |z^r^ A ) 
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